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Abstract 

We derive the flow equations for the symmetry unbroken phase of a dilute 3-dimensional 
Bose gas. We point out that the flow equation for the interaction contains parts which are 
non-analytic at the origin of the frequency-momentum space. We examine the way this 
non-analyticity affects the fixed point of the system of the flow equations and shifts the 
value of the critical exponent for the correlation length closer to the experimental result in 
comparison with previous work where the non-analyticity was neglected. Finally, we em- 
phasize the purely thermal nature of this non-analytic behaviour comparing our approach 
to a previous work where non-analyticity was studied in the context of renormalization 
at zero temperature. 
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1 Introduction 



Thermal effective actions are in general non-local in coordinate space because the temperature- 
dependent Green functions contain parts which are non-analytic at the origin of the momentum- 
frequency space [|1],@]. In a theory with two interacting scalar fields, one integrates out one of 
them to find the effective action for the other. At finite temperature, provided the coupling 
is weak, one usually proceeds by applying perturbation theory, and then making an expansion 
in powers of frequency and momentum in order to obtain a local effective Lagrangian. It is 
this latter expansion which leads to results which are not uniquely defined but depend on the 
path on the frequency-momentum plane through which the origin is approached. For example, 
when the perturbation is truncated at the self-energy level, the self-energy is non-analytic at 
the origin. The reason is that the expansion is around a singularity || . 

This effect was noticed for the first time by Abrahams and Tsuneto in the 60's, in the 
context of BCS theory, while they were studying time-dependent Ginzburg-Landau theory near 
zero temperature and near the critical temperature Q. Later it became clear that it is the 
origin of Debye screening and of plasma oscillations in QED These two different physical 

phenomena correspond to two different ways of approaching the origin of the momentum- 
frequency plane. The effects of the non-analyticity have also been studied in QED 3 and 
in QCD p|-|T3||. The non-analyticity is also present in the graviton self-energy ||14|, p"5|j and 
in higher-order graviton diagrams |l(J. Even in the much simpler case of interacting scalars 



the non-analyticity of the self-energy persists ||IT|-|2"0|| . In the case of interacting scalars, an 
interesting remark is that, whenever the internal propagators in a loop have different masses, 
the self-energy is analytic at the origin ||21|| . The reason is that the singularity is not at the 
origin anymore, allowing thus a uniquely defined expansion around the origin. 

This paper is based on the simple observation that an essential step of the renormalization 
group method (RG) applied to a theory with a self-interacting field, is to split the field in slow 
and fast components, and integrate out the fast field obtaining thus an effective action for the 
slow field |22|,[23[]. Therefore, according to the above discussion about thermal effective actions, 
when RG techniques are applied in the context of thermal field theory, we anticipate that effects 
originating in the non-analyticity will arise. 

We choose to examine this aspect of RG in the context of a 3-dimensional homogeneous 
self-interacting bosonic gas with weak repulsive interactions and discuss its possible physical 
significance in this case. However our analysis and conclusions should hold whenever RG is 
used at finite temperature. This choice of system was motivated by the renewed interest in the 
Bose-Einstein condensation due to its recent experimental realization. For the interacting gas, 
the approach which is most often used is that of Bogoliubov. However, this is just a mean-field 
type method and, in principle, one can improve upon it by using more sophisticated techniques. 
One possibility near the critical region is the renormalization group p4|-27 . 

In the case of the homogeneous gas, there is an extra, more important reason for looking for 
alternatives to the Bogoliubov approach. In the critical region, the Bogoliubov theory simply 
does not work because there are fluctuations around the mean-field that cannot be treated per- 
turbatively. This happens because, as the temperature approaches the critical temperature T c , 
the thermal cloud density develops an infrared singularity and thus diverges as the momentum 
tends to zero |28|. f2~9| . 

In section 2, we introduce the basics of the BEC formalism above the critical region. We 
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then apply Wilsonian renormalization and derive the flow equations for the parameters of the 
Lagrangian. We point out the non-analytic structure of the RG correction to the interaction 
term (vertex) and follow this non-analyticity as it propagates to the flow equation for the 
interaction. 

In section 3, we calculate the non-trivial fixed pont of the system of the flow equations and 
find the critical exponent for the correlation length. We examine the way the non-analyticity 
affects the fixed point. We note that taking the non-analyticity into account shifts the value of 
the critical exponent closer to the experimental result. 

In section 4, we compare our work with |3(| where the issue of non-analyticity in the context 
of renormalization is also discussed. We point out that the conclusions of j30]| hold only at T = 
whereas the non-analytic behaviour which we are investigating in this paper is purely thermal 
and vanishes at zero temperature, thus being completely independent of the non-analyticity 
discussed in . 

In section 5, we present our conclusions. 



2 Non-analyticity in the uncondensed phase 



When the two-body collisions between bosons are taken to be low-momentum or s-wave, the 
path integral representation of the partition function of the homogeneous interacting Bose-gas 
is given by 



Z(ji, (3, V, g) = Tre-^ 6 -^ = J 



-S[<j>,4>*] 



with the action 

S [0,0*] 



~ / dr [ d D x 

a Jo JV 



d 



h 2 



l 



(r,x)[^ - ^V 2 -/#(r,x) + -s|<Kr,x)| 4 



(2) 



In the low-momentum approximation the interparticle interaction can be described by the 
zero-momentum component of the Fourier transform of the two-body interaction potential. 
Thus, within this approximation, a repulsive, short-range potential can be characterized by a 
positive interaction strength g. In three spatial dimensions this interaction strength is related 
to the positive scattering length a of the interparticle interaction by the familiar relation g = 
Anh 2 a/m. The chemical potential is denoted by fi. The case /i < holds for T > T c and 
corresponds to the uncondensed phase whereas fi > describes the condensate which is formed 
when T < T c In this paper we will deal only with the uncondensed phase. Starting from (Q) 
we can derive the renormalization group equations for g and fi. This set of coupled differential 
equations can then be used for the study of universal as well as non-universal properties of the 
gas p4^p7|. In the following we will set h — 1. 

In order to implement the first step of the RG procedure (Kadanoff transformation), we split 
the field 4>(x) into a long- wavelength component 4><{x) and into a short- wavelength component 
&/>>(x). The short- wavelength field involves Fourier components which are contained only in an 
infinitesimally thin shell in momentum space of thickness A(l — dl) < |p| < A near the cutoff 
A, whereas the long-wavelength field has all its Fourier components in the sphere whose center 
is at the origin of the momentum space and its radius is A(l — dl). We impose no cutoff on the 
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frequency and apply the Wilsonian technique of consecutive infinitesimal shell integration only 
to the momentum and not to the frequency |]32| . 



We denote the volume of the shell by 5V P , the volume of the sphere by V p . The coordinate 
space volume is denoted by V. For simplicity we will be referring to as the lower or slow 

field and to 5<p > (x) as the upper or fast field. Whenever more compact notation is required we 
will be making use of the following: 

x=(t,x), p= (pq, p) with Pq = 2vm//5, 

\ dx = t dT Iv rf3x ' / dp = jJL I W? ' 

We integrate out the upper field and are left with an effective action for the lower field 

S eff [0<,0<] = ^[0<,0<] + ^TrLn[l - G>t\. (3) 



with 



For details on the derivation of this result and the approximations involved see The Tr 
denotes the trace in both the functional and the internal space of G > S whereas the tr denotes 
the trace only in the internal space (see below). The hat denotes that the corresponding 
quantity is a Schwinger-Fock operator |33|], 

U [P) { B*{p) J 

and 

r(rt = 9 - ( ^(x)^) 20<(x)0<(x) \ 

2\ 20* < (x)0* < (x) 40* < (x)0<(x) ) 

1 " ^ 

B(p) = B(p ,p) = - E (p) = ?--^ (4) 

tpo + -c/(p) 2m 

Note that the expression for S contains the coupling g. This enables us to perform a pertur- 
bative expansion over g in (§) in order to calculate it explicitly. We truncate this expansion to 
second order in g 

TrLn[l - G > ±] w Tr[-G>S - -(G>E) 2 ]. 

The first trace is: 

Tr[G>S] =fdx j dptv[G > {p)H(x)] 

= fdx \M^9 / + 2N[E(p)]] (5) 

J JSVp [In) 6 

where N[E(p)] = [e^ E ^" — is the Bose-Einstein distribution. We note that the first trace 
is quadratic in the modulus of the lower field and can therefore be interpreted as a correction 
to the chemical potential 
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df*=-gf -£^[l + 2N[E(p))). (6) 

The second trace is: 

Tr[G > £G > £] = J dp J dk J dx J dy e i{p - k){y ' x) tr[G > (p)S(x)G > (A;)S(y)] = 



dp J dk J dx J dy e i( - J 



i{p-k){y-x) 



4 

[16 B{p) B{k) </>%(x) 0<(x) ^MMv) + 4 B{p) B*{k) 0<(x) 0<(x) ^{y) + 
4 £T(p) B(k) <f>%(x) 0* (*) </><(</) 0<(y) + 16 B*(p) B*(k) ^(x) 0<(z) ^(y) Mv) } ■ 

(7) 

In order to simplify the above expression we change variables as follows: 

1. In the second and fourth terms in the square brackets, p — > — p and k —>■ —k. 

2. In the second term, x — > y and y — *■ x. 
The second trace now becomes: 

Ti[G > £G > £] = J dp J dk J dx J dy j 

| [ e i(p-k)(y-x) + e i(k- P )(y- X ^ B ^ lg ^* ^ 

+2e <<p-*)<*-.) 4 * <(x) * <(x) 0<(y) 0<(y) | (g) 

Changing variables again, k — > g = — p, yields: 



Tr[G > SG > E] = J dq J dx J dy ^ 

{2 Ji( g ) 4 0* (x) 0<(y) 0<(y) 



+ 

where 



e -iq{y-x) _|_ e iq(y-x) 



J 2 (q) 16 0<(x) </>*(y) 0<(y) } (9) 



•W, q) = / dp 5*(p) S(p + g) and J 2 (g m , q) = | dp £(p) B(p + q). (10) 

We note that p is the momentum of the upper field (integrated over the infinitesimal shell 
around the cutoff) whereas q is the momentum of the lower field, Fig.|l|. 

It is essential in the RG procedure, and in particular in the Kadanoff transformation, to 
recast the effective action obtained after integrating out the upper field in the form of the 
original action ([|). The first trace is in a form that can be interpreted as a correction to the 
chemical potential. This is not the case however for the second trace; there are quartic products 
of fields but, unlike the four-field coupling term in the original action (FJ), these are non-local in 
coordinate space, thus not allowing the effective action to be recast in the form of the original 
action. In other words, though we start from an action containing interactions which are local 
in coordinate space, the RG procedure generates more general, non-local interactions. 
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This is a well-known feature of RG, namely to generate extra terms that do not appear in the 
original action and have a more general form in comparison to what we started with |22],|3j].|35 



Provided that these extra terms are not relevant they can, in most cases, for the purpose of 
calculating universal properties, be discarded. 

In our case we can Taylor expand <p<{y) around <p < (x). If we truncate this expansion at 
leading order, 0<(y) ~ <f) < (x), we remain within the family of local interactions we started with. 
This procedure is usually called derivative expansion and is physically relevant only when the 
lower field is slowly varying both in space and in time. 

Wilsonian renormalization is compatible with the derivative expansion. The reason is that 
in Wilsonian renormalization we are interested in constructing an effective action for the slow 
field. This compatibility can also be seen from a more technical point of view; the derivative 
expansion of the lower field is equivalent to an expansion of Ji and J2 in powers of go an d 
q (this is easily seen from doing integration by parts). This means that truncating the 
derivative expansion at higher than the leading order would give momentum and frequency 
dependent corrections to the interaction. However such terms are known to be irrelevant (see 
for example ||34|| , page 128) and can therefore be omitted. 



The second trace becomes: 

,2 



Tt[G > ±G > ±] = J dx^j 

ilmpe^Jxiq)) (?)4>< (x)4><(x)4><(x) + lim[(e^ + e~^) J 2 (q)} 160* < (x)0 < (a;)^(x)0 < (x) 



'11} 



which is a local expression that can be interpreted as a correction to the coupling term of the 
original action 

dg = ~9 2 { m lim (0 MQo, q)l +4 hm [J 2 (g m ,q)]} . (12) 
[ (g£Mq|)->(o,o) (C>NH(°>°) J 

It is at this point that the non-analyticity enters our discussion. Because we are at finite 
temperature, the integrals over frequencies in J\ and J2, h and 1% respectively, become sums 
which can be easily calculated when we turn them into integrals on the complex plane through 
Poisson summation 

■W, q) = L 7S3 q> p) and -W. q) = / 7S3 W> q, p) (is) 



where 



h(% , q, P) - ^ B fa, p)B(p + g , p + q)) - E{p - q) - E{p) _ — , (14) 
h(q , q, P) - j ^B(p ,p)B(p + q , p + q)) - ^ (p + q) _ E{p) + iq rn ■ (15) 
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27rm//3. In the following, we will suppress the 



We have set expi/3g™ = 1 because q™ 
superscript of q™ for simplicity. 

The first sum, Ji, is non- vanishing at T = and is known as the regular term. The second 
sum, I 2 , is purely thermal and is usually called Landau term in the context of thermal field 
theory @. We observe that the successive limits of Ji(q ,q) coincide, i.e. 



lim lim Ji(g ,q) = / 



d 3 P 



?o^u| q |-*o - Jsv p (27r) 3 2E(p) 

whereas the successive limits of ^(SojQ) do not 



[l + 2JV[S(p)]]= lim lim Ji( 9o ,q) 
|q|->Ogo-+0 



lim lim J 2 (q , q) = ? lim lim J 2 (q , q) = / -^P ^WjV 2 ^)]. 

go-»0|q|->0 |q|-»Ogo->0 J5V P {2,7() 6 



(16) 



(17) 



The reason these two limits do not commute is that J 2 has a singularity at the origin 
of the momentum- frequency space ||. Of course, in the evaluation of the above limits, we 
interchanged the limits with both the integration over the momentum |p| and with the angular 
integration over 9, so our conclusion is not entirely reliable so far. In principle, one should 
perform the integrations over |p| and 9 first, and then take the limit. Unfortunately, in our 
case, the integration over |p| cannot be done analytically. However, we can perform the angular 
integration over 6 analytically before evaluating the limits, provided that we split the integral 
as follows 



d 3 p 



N[E(p) 



sv P (2tt)3 E(p + q) - E(p) 



«9o 



d 3 



p 



N[E(p + q)} 



5V P (27r)3£(p + q)-£(p) 



no 



(18) 



and perform the change of variables p — > — p — q in the second term, eliminating thus the 
dependence of the Bose-Einstein distribution on the angle 9 and making the angular integration 
possible. This procedure yields the result 



■M?o,q) 



d a 



P 



<5Vp (27r) 



+ 



1 



E(p + q) - E(p) + iq E(p + q) - E(p 



N[E(p)}. (19) 



It is crucial to note that this change of variables is not permissible in case we interchange the 
limit lini| q |^ hm 90 ^ with the integrations, because this causes both the terms in ( |T8"D to diverge. 
These two divergencies canceled each other before the change of variables p — > —p - q 0. 
Keeping this remark in mind we now perform the angular integration and find 



m 
4^ 



A-dA |q| 



(E + - E(p)f + ql 



(20) 



_{E_-E{p)f + ql_ 

where E + = E(\p\ + |q|) and E_ = E(\p\ — |q|). Instead of just taking the two successive limits 
in J 2 as we did before, which in the momentum-frequency plane corresponds to approaching 
the origin in the direction of the one or the other axis, we could approach the origin through 
any other curve, for example in the direction of any straight line go — a|q| . Here, of course, we 
should not forget that the frequency is discrete whereas the momentum is continuous. However, 
for the purpose of better illuminating the structure of J 2 around the origin, we shall make the 
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approximation that the frequency is continuous so that go 
Applying this parameterization to ([ 



lim J 2 (a\q\, |q|) 

|q|->0 



and then taking the limit |q 
d\p\N[E{p 



a|q| can hold for any real a. 
yields 



m 



2 


P 


2 


IP 


2 . 


h m 2 a 2 



(21) 



4tT 2 JA-dA 

which reproduces the first limit of flT7|) for a — > oo. This result was derived from ( p0|) and 
therefore is also not valid when the limit a — > is interchanged with the integration over |p|. 
However, if we do an integration by parts, we find 



lim Jgfalql 

|q|-0 



2tt 2 



A-dA 



,N) = 
rf|p||p| 2 



m 
2^ 



N[E(p)} 





P 


3 




2 


+ m 2 a 2 



A-cZA 





P 


2 


lp 


2 


+ m 2 a 2 



/Je^iV^p)] 



2m 3 a 2 



|p|2 _|_ m 2 fl 2j2 



tf[J5(p)] 



(22) 



We note that the surface term vanishes, as it is evaluated at the cutoff. This result not only 
reproduces the first limit of ( |TTD for a —>■ oo but also agrees with the second limit of flI7|) for 
a — > 0. If we perform the angular integration and apply the same parameterization to J\, at 



the limit |q| — > 0, J\ is independent of a and given by (0). Expressions fllED and ( p2[ ) are to 
be substituted in the correction for the coupling constant ([12]) . 

Before we proceed to the second step of the RG formalism, we parameterize the momentum 
according to |p(7)| — Ae~'. The purpose this change of variables serves is simply to make the 
flow equations more elegant. 

So far, the flow equation for the chemical potential is: 



\3 p -3l 

-g I ±^{l + 2N[e A e- 2l - f ,}} 



2tt 2 



(23) 



and the flow equation for the coupling is: 

r A3 e -3« 1 



dg 
dl 



-17 = ~9 



Ah 



-31 



2vr 2 



2tt 2 2[e A e- 2l -i4 
A 2 e- 21 



A 2 e 21 + m?a 2 



1 + 2N[e A e~ 21 - //]]+ 
(3N[e A e- 21 - //][! + N[e A e~ 21 - //]] 



2m 3 a 2 



[A 2 e~ 21 + m 2 a 2 ] 2 



N[e A e- 21 - fi] 



(24) 



where e A = A 2 /2m. 

At this point we apply the second step of the RG procedure, namely the trivial rescaling, 
whose purpose is to bring the effective action in the form of the original one. There are two 



q 



q\e l in order to re-establish 



stages, first, we rescale the momentum according to |q| 
the original cutoff A. Then we require that the effective Lagrangian has the same form as 
the original Lagrangian. This induces the trivial rescaling of the parameters of the effective 
Lagrangian. 



V -> V(l) = V e" 3 \ 
13 - (3{l) = (3 e- 2 \ 
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- 0(0 = 0e 3 '/ 2 , 

H -> /x(Z) = + A/i) e 2 ', 

5 -> #(0 = (<? + &g) e~ l . 



(25) 



The trivial rescaling of (3 implies that the frequency is rescaled as go — > 3o(0 e anc ^ therefore 



a = 9o/|q| a(l) e" 



(26) 



Recasting (123) and ( [Zip in terms of rescaled variables yields the flow equations for the corre- 
sponding running quantities 



and 

dg(l) 
dl 

where Ni 



-9(l)-9 2 (l) 



A 3 



dfj,(l) 
dl 



1 + 27V, 



2»(l)-g(l)^[l + 2N l 



2vr 2 \2[e A -fx(l)) 



+ 4 



A 2 



A 2 + m 2 a 2 (/) 



/3(0M[l + iV,] 



2m 3 a 2 (0 



[A 2 + m 2 a 2 (/)] : 



(27) 



3 /3(0[e A - M (0]. 



(28) 

1] 1 is the Bose-Einstein distribution in terms of the rescaled variables. 



3 Fixed Point 

It is important to investigate whether the path-dependence of the flow equation for the coupling 
has any consequences on quantities of physical interest. 

We look at a universal property, the critical exponent for the correlation length. This is 
calculated from the coupled system of (|27|) and fl28[) . In fact, in order to have an autonomous 
system, we should also take into account the flow equation for the inverse temperature, 

d -f = -2/3(0 (29) 



which is just a differential expression of the trivial scaling of (3 (25) . We observe that, although 
j3 appears in the equations for \i and g, these do not couple back to the equation for (3. We also 
note that the fixed point for (3 is zero, % = 0. The fact that /3* — complicates things because 
(3(1) appears in the flow equations not only explicitly but also through A^. This means that, 
when we evaluate the fixed point for the system of the flow equations, the right-hand side of 
(|27| ) and ( p8|) will diverge, because Ni diverges for /?* = 0. This problem is circumvented when 
we define a scaled running chemical potential M(l) and a scaled running coupling constant G(l) 
such that the set of equations for these new parameters decouples from the equation for j3\ 

M{1) = (3 A fi{l) and G{1) = A 3 (3 A g{l)/b{l) (30) 

where (3 A = m/A 2 , e> = 1/2 and b(l) = (3(1)/ (3 A is the scaled inverse temperature. In terms of 
these new, dimensionless parameters 



9 



dM(l) 

dl 
dG{l) 

dl ' 
db(l) 
~dT ~ 



-- 2M{1) 
G{1) 



2tt 



;G(l)b{l)[l + 2N l ] 



2tt 2 



G 2 {l)b{l) 



l + 2iV ; 
2[e>-M(0] 



+ 4 



A 2 



A 2 + m 2 a 2 (/) 



&(Z)iVi[l + iVi] 



2m 3 a 2 (l) 



(3 K [A 2 + m 2 a 2 ( 



-2b(l). 



(31) 



In the neighbourhood of the fixed point, the rescaled temperature is high and the approximation 
Ni « [(3(1) [e A - MO]] -1 = [6(0 [ e > - hoids pi. This yields the equations 



dM(l) 

dl 
dG{l) 

dl ' 



2M(l) 



2tt 



1 G(l) 6(0 [1 + 2 



6(0 



+4 



2tt 2 
A 2 



6(0 [c 



+ 



M(or 
i 



d6(Q 



A 2 + m 2 a 2 (0 [e>-M(0] 2 
-26(0- 



2[e>-M(0] [e>-M(0] 2 
1 2A 2 m 2 a 2 (Q 
[A 2 + m 2 a 2 (/)] 2 



+ 4 



A 2 



6(0 



A 2 + m 2 a 2 (0 e>-M(0 
1 1 



M(0 



(32) 



At the fixed point (M*, 6* = 0) the left hand side of (|32| ) is zero by definition. The form of 
the right hand side depends subtly on whether a is zero or not as we will see. Near the fixed 
point, the second term is dominant in the square brackets of the flow equation for the chemical 
potential in fl52]), 

' -2M{l)-—G{l) — -. (33) 



dl "" w 7r 2 " vv e> -M(0' 
For a = 0, we recall (p6[), which means that a(/) = for any /. Consequently, the equation 



for the coupling reduces to 
dG{l) 1 ~. 



c(0 



2tt 2 



G\l) 



6(0 



M(0] 



+ 



1 



M(l)f 



+ 4- 



6(0 



M(0 



+ 4. 



M(/)] 5 



(34) 

Recalling the trivial scaling of 6 (25), we see that the second and fourth terms in the curly 
brackets above are dominant near the fixed point, 



dG(l) 
dl 



G(l) ~ ^G\l) 



2vr 2 



M(IW 



(35) 



and calculate the non-trivial fixed point 



We linearize the system of fl3"2Q around the fixed point and find the largest eigenvalue A + = 
1.878. Therefore the critical exponent for the correlation length is v — 0.532, which agrees with 



the one found in 24 
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For a ^ 0, according to the trivial scaling of a (PB|), a* = oo and therefore the first and the 
third terms in the curly brackets of (|32|) vanish near the fixed point, as in the case of a = 0. 
From the remaining terms the fifth vanishes and, significantly, the fourth is also vanishing near 
the fixed point, leaving as dominant contribution only the second term. Consequently the flow 
equation for the coupling reduces to 

^-^-^K^mW (36) 

and the non-trivial fixed point is 

(M*,aA) = (^y,0). 

Linearizing around the fixed point we find that A + = 1.561 and therefore v = 0.640. It is 
interesting to note that we could have found the same fixed point and critical exponent directly 
from (|32|) , had we set a = oo and therefore a(l) = oo for any I. 

This situation is similar to what happens, for example, in the case of thermal QED for 
the photon propagator. Because the photon self-energy is non-analytic at the origin, different 
ways of approaching the origin lead to different dispersion relations and give rise to different 
types of excitations |T|, @. For short wavelengths, the dispersion relation is g 2 = q 2 + m 2 p , 
where m 2 p = (e 2 /2)(T 2 /3 + /i 2 /7r 2 ) is the thermal mass for the transverse photons whereas 
the longitudinal photons do not propagate. However, for long wavelengths, the transverse 
photons have the dispersion g 2 = up + |q 2 and the longitudinal photons have the dispersion 
q 2 = up + |q 2 , where u 2 p = |m|, is the plasma frequency at order e 2 . The phenomenon 
which we are describing here is of the same mathematical nature, the difference being that it 
is occurring not in the propagator but in the vertex between four bosons. To be more precise 
it is the vertex graph corresponding to J2, Fig.|l|, that exhibits the same singular behaviour as 
the photon self-energy in QED. 

Comparing the two critical exponents derived above with the known experimental value 
v = 0.670 [j36H , it becomes clear that the second procedure gives an improved estimate. In fact, 
our estimate is better even than the value v = 0.613 calculated in |24]] with the inclusion of the 
marginal three-body scattering term in the action. 



Flow equations and non-analyticity at zero and finite 
temperature 



Non-analyticity in the context of RG has been discussed before by Shankar in ||30|1 . In this 
work, the author gives a detailed overview of the RG approach to interacting, non-relativistic 
fermions in one, two and three dimensions and in certain instances (pages 161, 166, 170, 178) 
refers to the non-analyticity (or lack thereof) which appears in the one-loop RG corrections to 
the quartic interaction among fermions. This is highly reminiscent of the case we are studying, 
the essential difference being that we are dealing with bosons instead of fermions. This would 
render the main point of this paper - the study of a non-analyticity in the flow equation for the 
interaction - rather trivial and expected by extending |30|] to bosons. 
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This is not the case however, the non-analyticity we are studying is of completely different 
nature from the one studied in ||30|| . The RG calculations in |30| are at T = whereas ours at 
and the non-analyticity we are referring to is essentially thermal and vanishes at T = 0. 
To further clarify this point, let us consider the "zero sound" (ZS) graph which is studied in 



Eq.(315), page 161 of g(i 



4.1 Zero Temperature Non-Analyticity 
4.1.1 Zero Sound Integral 

In the zero-sound calculation of |j30 |, the following integral appears: 



La 2vr Loo 2vr [iu - k][i(u + Q) - k - q] (3?) 

where Q, q are the external frequency and momentum and u, k are the internal frequency 
and momentum respectively. The external momentum is constrained by the same cutoff as the 
internal momentum, —A < q < A. 

We focus on the integral over uj. When the external momentum and frequency are zero, the 
integrand has a double pole and therefore Si[0, 0] = 0, Fig.^. 

For non-zero external frequency and momentum, the integrand has two single poles. Let us 
assume that k > 0. If k + q > 0, both poles are in the lower half-plane and closing the contour 
from above yields Si[f2,g] = 0. However, if k + q < 0, the two single poles are in different 
half-planes and closing the contour either from above or below yields 



A dk 



—A 2ti n + iq 



Assuming k < we can argue the same way, so, only when k and k + q have different signs, 
the cu-integral is non-zero, Fig.|3[ 

Because of the fc-integration, there is always a range of k values for which k and k + q have 
different signs. Therefore any non-zero q results in a non-vanishing 



/A A lr n 
-A 2tt i I + IQ 



r A dk i 

iq 

which depends sensitively on how the limit {Q, q} — > {0, 0} is taken 



lim lim Si[Q, q] = ^ lim lim SAQ, q] = —1. 
4.1.2 Renormalization Integral 

In the RG calculation of |3(J| , a one- loop correction to the interaction is of the form: 

dk r°° du 1 



s 2 [n,q] 



-A+dA rA 
+ 

A JA-dA 



2-7T J-oo 2n [iuj — k] [i(u + Q) — k — q] 
As before, for zero external frequency and momentum, 5*2 [0, 0] = 0. 



(38) 
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For non-zero external frequency and momentum, we note that S2 differs from Si only in the 
range of integration of the momenta, k G [—A, — A + dA] U [A — dA, A] and q G [—A + dA, A — dA] . 
This means that \k\ > \q\ and consequently k,k + q have always the same sign throughout the 
integration over k. Therefore S2IP, <?] = which has no dependence on how the limits of the 



external frequency and momentum are taken [30|. 



The conclusion is that, when performing RG calculations, non-analyticities at the origin of 
external frequency and momentum space vanish, even when they are present in the correspond- 
ing zero-sound calculations. 

4.2 Thermal Non-Analyticity 
4.2.1 Zero Sound Integral 

Now consider what happens at T ^ 0. In zero-sound calculations one has to perform the 
integral 

SliV m ,q] = [ A f-\ if ^ (39) 

J -A 2n (3 n fz' 0O [iuj n - k\ [t(un + il m ) — k — q\ 

where u> n = (2n+l)7r//3, fl m = (2m+l)7r//3 are the discretized internal and external frequencies 
respectively 

When the external frequency and momentum are zero, we find 

s 'i - ' = -/ A A Sr osh "ifi-- L < 40 > 

For non-zero external frequency and momentum, we obtain 

which gives the correct zero temperature limit Sj^°[Q m ,q] = Si[Q,q]. At T 7^ 0, 

lim lim S7\n m , q] = ^=3 0, (42) 



lim lim SR^g] = - f\ ^cosh" 2 ^] ^° -1 (43) 

and therefore the non-analyticity of the zero-temperature zero-sound calculation persists at 
finite temperatures, Fig.|j. 

4.2.2 Renormalization Integral 

As at T = 0, the integral appearing in RG calculations, differs from Sj only in the range over 
which the internal and external momentum are integrated. 

When the external frequency and momentum are zero, we find 



dk (3 __ 2 fc/? T ^o 



S 2 T [0,0] 



-A+dA pA 
+ 

A JA-dA 



^cosh- 2 [M]^ . (44) 
2vr 4 1 2 J v ' 
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For non-zero external frequency and momentum, we obtain 
Sj[tt m ,q] 



r-A+dA 


f 1 


L + 


JA-dA 



dkl ~ 



iu n - k] [i(u n + Q m ) -k-q}' 



(45) 



Poisson summation yields, Fig.[|. 
Sj[tt m ,q] 



p-A+dA 


/ 1 


L + 


JA-dA 



dk 



2ir Q m + iq 2 [ 
which gives the correct zero-temperature limit 

) [Q m ,q] = 0. 



tanh 



.k/3. 



tanh 



.(k + q)f3 



(46) 



However, unlike the zero temperature renormalization integral, the finite temperature renor- 
malization integral is non-analytic, because 



lim \imSj[n m ,q] = 0^0, 



lim lim Sj{tt m ,q] 



-A+dA 



+ 



A 

A-dA 



^^cosh- 2 [M]I^ . 
2vr 4 L 2 J 



(47) 
(48) 



Therefore, at finite temperature, the non-analyticity of the zero-sound one-loop graph persists 
but there is also an extra, purely thermal non-analyticity appearing in the finite temperature 
renormalization one-loop graph. 

The zero temperature non-analyticities are due to the splitting of a double pole into two 
single poles residing in different half-planes. The finite temperature non-analyticities are due 
only to the splitting of a double pole and appear even if the resulting two single poles reside in 
the same half-plane. Therefore it is only natural that extra non-analyticities appear at finite 
temperature in addition to those existing at zero-temperature. 



5 Conclusions 

We have investigated thermal non-analyticities at the origin of the momentum-frequency space 
in the context of Wilsonian renormalization. We have shown that taking them into account leads 
to an improved estimate for the critical exponent of the correlation length, when approaching 
the critical region from the symmetry unbroken phase. It is well-known that the e-expansion 
gives values closer to the experiment than ours, however one should remember that the e- 
expansion is asymptotic . It is therefore crucial to find safer ways of calculating the critical 
exponent, such as the improved momentum-shell method which was used here. 

When approaching the critical region from the symmetry-broken phase there is an excellent 



estimate v = 0.685 [24]. It may be possible to improve this value by taking into account the 
non-analyticity as we did for the approach from the symmetry unbroken phase. Work on this 
issue is currently under progress. 

Finally, we have pointed out that when one applies Wilsonian renormalization at finite 
temperature, one may encounter non-analytic behaviour which is different in nature from the 



non-analyticity (or lack thereof) discussed in |30 
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p 





Figure 1: The two contributions of the RG correction to the interaction g, J x and p is the 
frequency-momentum of the upper field (the momentum p is integrated over the infinitesimal 
shell around the cutoff), q is the frequency-momentum of the lower field. 
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Figure 2: The integrand of the cu-integration has one double pole and the integral is zero. 
Any non-zero external momentum splits the double pole into two single poles. There is always 
a range of k for which the two single poles reside in different half-planes and consequently the 
integral is non-vanishing. 
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Figure 3: The integrand of the cu-integration has one double pole and the integral is zero. 
The introduction of non-zero external momentum splits the double pole into two single poles. 
However, because the fc-integration is over the infinitesimal shell near the cutoff and the external 
momentum q takes value below the infinitesimal shell, the two single poles are always on the 
same half-plane and the integral remains zero. 



20 





Figure 4: The integrand of the cu-integration has one double pole. The introduction of non-zero 
external momentum splits the double pole into two single poles and the integral is non-vanishing 
regardless of whether the two single poles are on different half-planes or not. 
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